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MULTI-STEP RUNGE-KUTTA METHODS 


SUMMARY 


The multi-step methods presented in this report differ from the classic, 
single-step Runge-Kutta process by utilizing the numerical results of the 
previous integrations. As developed here, the multi-step process avoids the 
complex Taylor expansions of the functions involved and instead employs the 
writer's previously developed method of quadratures. 


INTRODUCTION 


The classic Runge-Kutta method, which is a single-step process, has a 
number of pleasing properties, but since it does not utilize previous numerical 
results of the integration, its efficiency is impaired. This defect will be some- 
what remedied here by developing a multi-step method that is quite analogous 
to the single-step Runge-Kutta process. By utilizing the evaluations at pre- 
ceding points of the function in the differential equation to be integrated, some 
reduction in the number of evaluations required in the single-step process is 
expected. 

Perhaps the greatest obstacle in the development of a multi-step process 
is the tedious — and formidable — task of making Taylor expansions of the func- 
tions (those for both the current and preceding points) about the current point. 

In a recent article [ 1] , the writer has developed a "method of quadratures" 
for obtaining the system of equations that determine the parameters of the 
(single-step) classic Runge-Kutta process; this method obviates the need for 
making the complex Taylor expansions. In this report we will show how the 
method of quadratures can be extended to the multi-step process. 


THE TWO-STEP FOURTH-ORDER SOLUTION 


For simplicity, we will consider only the initial value problem 

y’ = 


f(x,y) y(x„) = y 0 


( 1 ) 



though the results may be extended to systems of differential equations. We 
define the sequence of equations 


ki (0) = hf (x 0 , y 0 ) 


= M (x 0 + a 2 h, y 0 + j3 21 k/°* 


( 0 ) - 


(■ 

^x 0 + a z h, 


Yo + 03 i k i 


(°) + 


^ 32 k 2 


( 0 ) 


(2a) 

(2b) 

(2c) 


where x 0 and y 0 designate the current point. 
Let 

x_! = x 0 - h 

y-i = y(xo-h) 


(3) 


refer to the preceding point separated by the step-size h. If we replace x 0 , y 0 
in equation (2) by x_ t and y_j (and make the appropriate changes in the super- 
scripts on the k’s) , we get 


( = hf 

T 

>> 

a 

o 

£ 

(4a) 

(-1) = hf 

^x 0 + (a 2 - 1) h, y_j + /? 2 i k i ( 1) j 

(4b) 

™ = hf 

(x 0 + (a 3 - l)h, y-!+ /3 31 k/ 1) + /3 32 k 2 1 ^ 

(4c) 


Inequations (2) and (4), and subsequently, the superscripts will distinguish 
the current from the preceding point. 


We will show that 

y(x 0 + h) » y 0 + + a> 2 ^k 2 ^ + w 3 ^k 3 ^ 

V’V- 1 * hV-‘V‘ I) + “ 3< " 1>kj< ‘ 1> 


(5) 


with appropriate values for the weights and the parameters in equations (2) 
and (4) , will give a solution of fourth-order accuracy at (x 0 + h) . Since 
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the evaluations of the functions in equation (4) will already have been made, 1 
only the three evaluations in equation (2) will be required to find the solution 
at the next point (x 0 + h) . Because these evaluations are the most time- 
consuming part of the computation and since the classic Runge-Kutta theory 
required four evaluations for this order, 2 the advantage of this multi-step 
method is evident. 

Suppose now that equation (1) is reformulated as 

, x o +h / v 

y(x 0 +h) - y(x 0 ) = J ^ f(x, y(x)J dx (6) 


in which form it includes the initial conditions. Now let us replace the right 
member of equation (6) by the approximating quadrature sum 


u> 


( 0 ) 


( 0 ) 


i f (x 0 , y 0 ) + w 2 ' f ( Xl , yi ) 


( 0 ) , . , 

+ co 3 'f(x 2 , y 2 ) 


L (7) 

+ w/ ^ f(x_ 1 , y-j) + cu 2 < ^ f (x- 2 , y- 2 ) + co/ ^ f(x_ 3 , y- 3 ) 


where for convenience we have set 


Xi = x 0 + a 2 h 

x 2 = x 0 + o; 3 h (8) 


and 


x_i = x 0 - h 

x _2 = x 0 + (a 2 - l)h 

X— 3 = x 0 + (a 3 - i)h 


( 9 ) 


and where y^ will be used for y (x^ ) . We will now show how we may replace 
each term in equation (7) by its corresponding term in equation (5). 


Except initially. Thus, in common with all multi-step methods, it is not 
self-starting. 

2 This assumes that the errors in the two cases are comparable — which need not 
be the case. 
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As we have shown in [ i] , in order for equation (7) to be a quadrature 
approximation of equation (6) with fourth-order accuracy in h, the following 
relations must hold: 


h(co 


(-1) (-1) (-1),' (0) (0) (Oh 

l w 2 w 3 i w i ^2 w 3 ) 


*1. \ 

X— 2 

x j 
x — 3 


x J 

x 0 

J 

x i 


a i 


Uj, j=0, 1,2,3 


( 10 ) 


where u^ is given by 


r X ° +h i, (x.+ h,^/ 1 . 

Uj= f xdx= j+1 , J-0,1,2,3 


(11) 


The reason for showing the partitions in the matrices in equation (10) is evident. 

If we substitute equations (8), (9), and (11) into equation (1©) and 
equate coefficients of like powers of h, we get the matrix equation that may 
replace equation (10) 3 


CO 


("I), ,(-!>, ,(-!) 


co 


CO 3 


(0) (0) (0) 


CO. CO 


co 3 


= (1 ~ 


1 (a 2 -l) (0:2-1 ) 2 (a 2 - 1)' 

i_ _L“r_ i_)_ 2 _(_"r il! 

100 0 


v 


«2 

«3 


) 


a 0 


«3 


CUn 


«3 


( 12 ) 


3 See reference [ 1] for more details about this operation. 
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These give some of the equations in the multi-step Runge-Kutta process. 

Using the quadrature method, let us now examine the further conditions 
by which equation (7) may be replaced by equation (5). If we compare analogous 
terms in equations (5) and (7) we see that the differences are in the expressions 
that replace y^ in the functions f(x i5 y^) ; the substitutions for y^, shown in 

equations (2) and (4) , substantially represent the unique aspect of the Runge- 
Kutta process. These interchanges for the y i in f(x^, y j) might be permissible 

if the terms in equation (7) are replaced by expressions of an acceptable order, 
the order three in this case, not four, because of the factor h in both equations 
(5) and (7). To attain this precision certain conditions, which will give other 
Runge-Kutta equations, must be imposed. 

Thus, the first term in equation (7) that will be affected by substituting 
replacements for y^ is f(x l5 y t ) , which is to be replaced by 

f^o + Yo + /*2i k / ^ (13) 

as may be seen from equation (2b) # This may be done by making a quadrature 
approximation for y t that replaces it with (y 0 + /3 2 ik 1 (^) ) , as is shewn in 
equation (13); to attain the requisite precision indicated previously, the 
following conditions must hold f 1] : 

(h/?2l)x 0 0 = U 0 (O! 2 ) 

(h/3 21 )x 0 = uj(a 2 ) 

(h/3 21 ) x o 2 = u 2 (a 2 ) 

(h/3 21 )x 0 3 = u 3 (q! 2 ) 

where 

x 0 + a i+1 h r 

u r ( a i+l> = J x x 
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Similarly, to replace f(x_ 2 , y_ 2 ) in equation (7) by f(x_ 2 ,y_ 1 + /3 21 k/ ^ ) , 

as shown in equation (4b) , with the necessary precision, we must also satisfy 
the conditions 

(h/ 3 2i ) (x„ - h)° = u 0 (o! 2 ) 

(¥21) (x 0 - h ) 1 = u^oia) 

(¥ 21 ) ( x o ~ h) 2 = u 2 ( 0 ! 2 ) 

(hjS 21 ) (x„- h) 3 = u 3 (a 3 ) (16) 

where 


^(“i+i) 


x o + (^i+i*" 1 ) 11 I x o + ( h l r + (xo-h) 1 ^ 1 

f , x r dx = ~r — — 

■'x 0 - h r + 1 


r - 0, 1, 2, 3 (17) 

As may easily be seen , equations (16) and (17) can be obtained simply 
by replacing x 0 by (x 0 - h) in equations (14) and (15); this is to be expected 
if we note that equation (4b) is obtained from equation (2b) by the same trans- 
formation. 

Similarly, we can replace f(x 2 , y 2 ) and f(x_ 3 , y_ 3 ) in equation (7) by 
the functions shown in equations (2c) and (4c) ; this is accomplished by re- 
placing y 2 and y_ 3 with appropriate quadratures. 

The conditions for replacing y 2 by (y 0 + /3 3l k/^ + /3 32 k 2 ^), as 
shown in equation (2c) and with the requisite degree of precision three, are 

(h0 31 )x o ° + (h^Xi 0 = u 0 (q! 3 ) 

(h0 31 )x o + (h/3 32 ) Xl = u^a-g) 

(h0 3 i)x„ 2 + (fy3 32 )x 1 2 - u 2 (a 3 ) 

(h/3 31 )x 0 3 + (h^g^xj 3 = u 3 (q! 3 ) (18) 

where u r (a: 3 ) is given by equation (15). 
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Similarly, to replace y_ 3 inf(x_ 3 , y_ 3 ) by (y_ft- 0 31 k/ lj + /3 32 k 2 ( **), 
as is indicated in equation (4c) , with the required degree of precision, 
we must satisfy 

(h/3 3 j) + (h/3 32 ) = u 0 (o; 3 ) 

(hiS 31 ) (x 0 - h) + (h/3 32 )x 0 = u 1 (q! 3 ) 

(h/? 3 i) (x 0 - h) 2 + (h/3 32 )x 0 2 = u 2 (Q! 3 ) 

(h /3 31 ) (x 0 - h) 3 + (h/3 32 )x 0 3 = u 3 (Q! 3 ) (19) 

where u r (Q! 3 ) is given by equation (17). As before, we can see that equation 

(19) can be obtained from equation (18) by replacing x 0 with (x 0 - h) — and, 
of course, xj with xq — with corresponding changes in the right members. 

THE RUNGE-KUTTA MATRIX EQUATION 


In [ 1] we have shown, for the classic single-step case, how the conditions 
we derived by using our quadrature method — in this present case the analogous 
equations (14), (16), (18), and (19) — are used to give an applicable Runge- 
Kutta matrix equation. 

We will here reconstruct the general idea without too detailed a treat- 
ment ( see also Appendix A) . Using the aforementioned equations, let us form 
the following weighted expressions: 

tof 1) u. (<* 2 ) + <^ 3 ( 1) W (a 3 ) + Co/ 0 * u. (a 2 ) + (" 3 ). j=0, 1,2,3 (20) 

Let us first dispose of the case j=0. We must first note, by equations 
(15) and (17), that 

“0 <“ 1 + i> * “ Vi h 

so that the first equations in (14) , (16) , (18) , and (19) are identically 
satisfied if 

a 2 = 021 

a 3 = 031 + 032 
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in general, the equations 


“i + i’'W + 'W + - • • + ' 3 i+i,i, 1=1 ’ 2 (N_1) (21) 

hold as they do in the single-step case. 

Let us illustrate the result of using equation (20) when j=2. Using 
equations (15) and (17) to replace the u. 's, we have 

(x 0 2 h) |cj 2 (_1) q^ + cu 3 ( 1) a 3 + u> 2 (0) Q! 2 + w^°^q! 8 J 

(2x 0 h 2 ) Jw 2 (-1) (-Q! 2 +-i-cv 2 2 ) + u) 3 (_1) (~a 3 +±a 3 2 ) + cu 2 (0) (|t* 2 2 ) 


+ aJ 3 0) (|*«3 2 ) J + ( h3 ) { w 2 ^ (j^ a 2 3 - a 2 + “2) 


+COo 


( -1) / 1 3 2 

! ra| - a, 2 


+ “3 ) + W 2 


(0) h„3 


“2 +CU 3 


( 0 ) 


'i-.s\ j 


( 22 ) 


The corresponding expression for (20) is as follows, 4 if we replace 
u 2 ( cn 2 ) , u 2 ( a; 3 ) and u 2 ( a 2 ) , u 2 ( a; 3 ) by the left members of the equations ( 14) , 
(16) , (18) , and (19) : 

(x 0 2 h) ju> 2 1q; 2 + w 3 _1 o! 3 + co 2 °a 2 + u 3 °a 3 J 

+ (2x 0 h 2 )J w 2 (_1) (-a 2 ) + w 3 (_1) (/l 32 a 2 -Q! 3 )+ u) 2 (0) . 0+ u> 3 (0) 03 S2 ce 2 ) j 
- ) 2~ 1 ^ a 2 + w 3 ^ (/ 3 32 a 2 2 - 2 Pz 2 0i 2 + a z ) + ^2^ • 0 


+ (h i ) 


+ w 3 (0) (P32 a 2 2 ) 


(23) 


If now we equate powers of h in equations (22) and (23) , we obtain the 
two Runge-Kutta equations: 

w 2 ^ a 2 2 + ^ ( a z' ” 2Pz2 a 2) + w 2 ^ a 2 2 + W Z ( a Z 2 ~ 2 Pz2 a 2) ~ ® 

(24) 


4 We omit the h term since it contributes nothing new. 
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<4 ^ (-« 2 S+ 3a 2 2 ) + <4 ^ + 3a 3 2 + 3/3 32 a 2 2 - 6(3 32 a> 2 ) ~ 4>°) o; 2 3 

+ 4 °> (3p 32 a 2 Z - Q!3 3 ) = 0 (25) 

No new results will be derived if we treat equation (20) similarly when 
j = 1, 3. 

Without repeating in detail what has been developed in [ 1] for the single- 
step classic Runge-Kutta process, we will note that we can replace (0) in 

equation (20) by c4^cr and, similarly, 4 ^ by 4 ^(at. - !)• A repetition 

of the previous analysis will yield another Runge-Kutta equation. 


We will see that, as in the single-step process, if we can find other sets 
of what we have called generalized Runge-Kutta weight coefficients to use in 
equation (20) we should eventually be able to produce all of the relevant Runge- 
Kutta equations. We will next show how we can get all of these coefficients 
from the single-step process and avoid the basic, but tedious, work in the 
illustration above. 


We will anticipate the results we would obtain in the form of a matrix 
equation: 5 


/ 


(- 1 ) 


CJ 3 


(-1) 


CO 2 


( 0 ) 


(-0 „ 
Cd 3 p 32 


0> 3 


( 0 ) 


\ 


/Wi\ ( -i) ( o) (o) 

td 2 ^ ( Q?2 - ^ ) w 3 ( (T 3 — 1 ) CO 2 a 2 w 3 «3 


( 0 ) _ 

<^3 P 32 


/-I 


<*2 


-al (j ce 2 - 1) 


\ 


P 32 a 2~ a \ CP 32 a 2 - o a 3 32 a 2 4 a l) 


1 2 


1 3 

■7 ®2 


y Ps2 a 2- 77 “3 (^32“2- 3 “3) 


(26) 


where 0, the right member of equation (26), is the appropriate null matrix. 


THE TWO-STEP FIFTH-ORDER RUNGE-KUTTA EQUATIONS 


The process for obtaining the generalized Runge-Kutta weight matrix, 
even for the relatively simple two-step fourth-order case, becomes somewhat 

5 As we will see later ,only four of these equations are relevant. The solution 
for the two-step fourth order is given in Appendix B. 
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involved when rigorously justified. We will now examine a procedure that is 
an extension of the one-step process, and that can easily be extended to any 
order, as is seen in [ i] . Let us demonstrate this by developing the two-step 
fifth-order Runge-Kutta matrix equation. With a slight modification of the 
results derived in [ 1] , we will first write out the one-step Runge-Kutta matrix 
equation for the fifth order with five evaluations. 6 


Let us define the matrices 



( 27 ) 


( 28 ) 


6 The exact number of evaluations is unimportant for our demonstration. 
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where 


n 




( 032 0 

0 

o\ 

n 

w 4 a 4 

n 

w 5 a 5 ) 

042 043 

0 

0 



^052 053 

£54 

•/ 

c<i> - 

0 




(j) 

c 3 

0: 32«2 J 




C, 0> - 

PiZ a 2 + 

043 ^3 



c 5 - 

052 a 2"* + 

053°^ + 054 

<y-l , 

j 


/ (n) (n) (n) (n) 
1 7 2 73 74 75 


n = 0, 1 


( 29 ) 


j = 1, 2, 3 


(30) 


and where 7 


,( 0 ) 


72 ’ (7j) ~ ( w 4^43 + £*32 + ( w 5^54) £*42 + 75^/3 52 

73 (0) (7 a ) = 74°^43 + 75^ 053 

74 (0) (7 i )= 7 5 (0) ^54 = 0 

7 5 ( ° ) ( 7 i )= 0 

The single-step fifth-order Runge-Kutta matrix equation with five 
evaluations 8 is given by 


n (t V°> = 0 


(31) 


7 The values for 


7- (0) (7 a ) 


(j = 2, 3, 4, 5) may be derived from y. 


( 0 ) 


J 


("j). 


defined by equation (29) , by replacing cn. in the latter terms with the 
corresponding y(^) . 

8 It has been shown [2, 3] that the fifth order requires at least six evaluations 
for a solution. Since we are using a two-step process, we would like to reduce 
the number of evaluations by one. 


11 



where 0 is the appropriate null matrix and where 


(0) _ r 

X “ C 


and C is given by equation (28) , 
get 


± 

2 

« 2 2 

1 

3 

«2 3 

1_ 

2 


3 

2 

«3 

3 

^3 

1 

9 

J_ 

„ 3 


3 

a 4 

_i_ 

o 

«5 2 

J_ 

3 

V 



( 32 ) 


Combining the matrices in equation (32), we 



/[ 

'iM [ 

-T “«'] [ 

4“^ 


\p3Z a 2 

[/ 3 32 Q! 2 2 

- y“3 3 ] [ft2“2 3 


x'°>= 

, [(^42 a 2 + 043 a 3 

*'| OJ 4 2 ] [ (^42 a 2 2 + 043 a 3 2 ) 

- y “4 3 ] [S42“2 3+ ^43 “3 3 



\ [(^52 a 2 + (^53 a 3 + @5i a i) 

-4-a 5 5 ] [ (/3 52 a 2 2 + ^53«3 2 + / 5 s 4 “ 4 2 ) - " 3 " “5 3 ] [(fe“2 3+ ^S3“3 3 + 

/ 3 54 “4 3 r^5 4 V 


(33) 


The two-step fifth-order Runge-Kutta process with five evaluations 
will take the linear form 


y(x 0 + h) 


5 

Yo + 2 


i=l 


«? k ( 


5 

+ 2 
i=i 


W <°> k ( °> 


(34) 


where 


k. (0) = hf^c 0 + a i h > Yo + 2 k j (0) j 


k[ ^ = hf 


i-i 

x 0 + (x -1) h, y_j+ 2 p k 

j=l J J 


(- i) 


(35) 


( 36 ) 


and where the superscripts (0) and (-1) refer to the current and preceding 
points, respectively. 
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To accommodate the antecedent point, certain modifications will have 
to be made in the two matrices in equation (31). In the case of the matrix 

Cl ^ , as we may see from equation (26) , this is quite simple because, aside 

from the appropriate superscript on o>., we must replace a. by (a. -1). With 

1 (- 1 ) 1 1 
these changes we will call the resulting matrix n . 

The accompanying change in the matrix in equation (33) is not as simple, 
and only the result will be given for the order of interest here. 9 Thus, if we 
write for the matrix in equation (33) 


(37) 

« 


(38) 

or, substituting the values for the elements from equation (33) and using the 
superscript (-1) to distinguish it from the matrix for the current point, we 
have 


9 This is developed in detail in Appendix A. 


11 

e 21 

e 3l 

'12 

e 22 

e 32 

*13 

e 23 

e 33 

*14 

e 24 

e 34 


we should replace this by 

f e li ( e 2 i “ 2e u ) (e 31 + 3e 41 - 3e 21 ) 

X = e i 2 ( e 22 “ 2e I2 ) (©32 + 3e 12 - 3e 22 ) 

1 e i 3 ( e 23 “ 2e 13 ) (e 33 + 3e 13 - 3e 23 ) 

\e 14 (e 24 - 2e 14 ) (e 34 + 3e 14 - 3e 24 ; 
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v (-l) 


032 a 2 


+ 043*3 ) 


+ a ?\ 

- J°3 j ^32“! 2 " J»3 3 - fe“2 + “3) 

-J“lj jj^«“2 , + P«“» -3“/) 

-2^32 ff 2 +Pi3 a 2 ~ 2 


<0S2“2 + Ps3 a 3 + 


0S4«<) -J“5 2 j 


(- 1 ^- !“*♦«*) 

- J“S + 3/3=2“2 
-f ■ “3 - 3^32“2 ! + “3 ! )j 


whose elements are partially shown in equation (39) . 


( 39 ) 


We can now combine our results to give us the matrix equation for the 
two-step fifth-order process. This can be written 


.(- 1 ) 


a 


(0) 


__x. 


(- 1 ) 

To") 


(40) 


where the partitioned matrices shown are as previously defined and 0 is an 
appropriate null matrix. 

As an illustration, let us return to the two-step fourth-order equations 
with three evaluations at each point, then the matrices in equation (40) become 


( 

-1) 


-1) 


,2 t°) 

„ (0) 

0>2 

W 3 


o> 2 

tu 3 

w 2 ( 

_1) (a 2 -l) 

W 3 ( 


-1) 1 

w 2 ( 0) a 2 


»2< 


0, 3 < 

' 1} (a 3 

-l) 2 f 

W2 <V 

-2 ,#) ^ 




0 


(0) a 
0 )j £32 

0 


*T 3 2( *2-1) 


0 


(0> 

W 3 032*2 

0 


^0 3 2( a 3~l) 


0 


«* ( °^32 a 3 

0 


0 


"‘W 

“ 2 - 1 ) 1 

0 

(Oh, 

032 


1 2 
-2*’ 


(p32 a 2~ 2 a i) 


-H) 


(032«2-J«3 2 ) 


(032®2 2 -f «3 3 

-Pn a 2 + a iJ 

(‘H 

^32°'2 2- 3 a 3 3 ^ 




\ 


r^32 tt 2 3 “ ^3) + 3 ^32«2" f *3*) 
-3032*2* + QT 3 3 J 


-±„j 


^22“2 S -J“i ij 


( 41 ) 
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which give both the Runge-Kutta equations for the fourth-order method 10 and 
the fifth-order terms from which we may find the truncation error. 11 If the 
latter terms are ignored, equation (41) will adequately represent the required 
equations. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 
Huntsville, Alabama, September 1, 1967 
039-00-24-0000 


i. 


10 These are the equations in which the sum of the exponents in the product of 
the a's and /3's in any term is always less than the order of the process. 

llr These equations, aside from some misprints, may be found in [ 4] . 
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APPENDIX A 


DERIVATION OF THE TWO-STEP FOURTH-ORDER EQUATIONS 


The conditions we have derived using the quadrature methods for the 
two-step fourth-order process are given by equations (14), (16), (18), and 
(19) . We will now display these equations in an expanded form in Table I and 
show how we can obtain from these the Runge-Kutta equations. 

In the first two groups of Table I are given the equations (14) and (18) 
for i = 0, 1, 2, 3, with the u (a 2 ) and u. (a 3 ) terms in the right members 

expanded by using equation (15). In the last two groups in Table I are the 
expanded equations (16) and (19) , where these are obtained from the first 
two groups by replacing x 0 by x 0 - h. 

In Table II are shown the results of combining similar terms of the left 
and right members of the equations in Table I. For convenience, in forming 
equation (20) we have associated each group in these tables with a weight 

coefficient o> ^ or ; and in Table I we have retained the origin of each 

equation by the appropriate label u. (a 2 ) or u. (o? 3 ). 

With the aid of Table II it is quite easy to form the expressions for 
equation (20) and to derive the matrix equation (26). 

It will be observed that by extending Table I the Runge-Kutta equations 
for any higher order multi-step process may be obtained in a similar manner. 
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TABLE I. THE CONDITIONS FOR THE TWO-STEP FOURTH-ORDER PROCESS 
GIVEN BY EQUATIONS (14), (16), (18), AND (19). 


¥21 = ha 2 

(0) (x o h)02i = (x 0 h)a 2 + -^h 2 0! 2 2 
2 

(x 0 2 h)|3 21 = (x 0 2 h) a 2 + (x 0 h 2 ) a 2 2 + — h 3 a 2 3 

3 ** j 

(x 0 3 h)jS 21 = (x 0 3 h)Q! 2 + (“^ x 0 2 h 2 ) a 2 2 + (xoh 3 )a 2 3 +-|-h 4 Q ! 2 4 


= Uo(a 2 ) 

= ^( 0 - 2 ) 
= u 2 (a 2 ) 
= u 3 (Q! 2 ) 


h(ft}i +P32) 

JO) x 0 hO 31 +/S 32 ) + h 2 (/3 32 u 2 ) 

x 0 2 h(j3 31 +/3 32 ) + 2x 0 h 2 (/3 32 a 2 ) + h 3 (/3 32 a 2 2 ) 

x 0 3 h(fti ¥ 32 ) + 3x 0 2 h 2 (/3 32 O! 2 ) + 3x 0 h 3 (/3 32 a 2 2 ) + h 4 (/3 32 Q! 2 3 ) 


= ha 3 

= Xoha 3 + -| h 2 a 3 2 

= x 0 2 ha 3 + x 0 h 2 Q' 3 2 + h 3 a 3 3 

= x 0 3 ha 3 + Y x 0 2 h 2 a 3 2 + x 0 h 3 a 3 3 +-i-h 4 a 3 4 


= u 0 (q: 3 ) 

= M« 2 ) 

= u 2 (a 3 ) 
= u 8 (q- 3 ) 


a 2 h = a 2 h 

(-1) (a 2 h) (x 0 -h) = (or 2 h) (x 0 - h) + h 2 a 2 2 

2 ^ 

( 0 ! 2 h) (x 0 - h) 2 = (a 2 h) (x 0 - h) 2 + (x 0 h 2 ) « 2 2 - h 3 a 2 2 + — h 3 a 2 3 

(a 2 h) (x 0 - h) 3 = or 2 h (x 0 - h) 3 + h 2 a 2 2 (x 0 2 - 2x 0 h+h 2 ) + h 3 Q! 2 3 (x 0 - h) + -“h 4 ^ 4 


= u 0 (a 2 ) 
= «>i(a 2 ) 
= Q 2 (« 2 ) 
= u 3 (a 2 ) 


a 3 h = a 3 h 

(a 3 h) (x 0 - h) + h 2 (jS 32 o; 2 ) = (ho; 3 ) (x 0 - h) + Jh 2 a 3 2 

(a 3 h) (x 0 - h) 2 + 2(x 0 - h)h 2 ((3 32 o; 2 ) + h 3 (/3 32 a 2 2 ) = (ha 3 ) (x 0 - h) ^h 2 ^ 2 (x 0 - h) + h 3 a 3 3 

(QJ 3 h) (x 0 -h) 3 + 3(xq - h) 2 h 2 (^ 32 Q! 2 ) + 3(x 0 - h)h 3 (/3 32 Q! 2 2 ) + h 4 (/3 32 a 2 3 ) = (ha 3 ) (x 0 -h) 3 


= Uo(“3) 
= Q l(«3) 

= u 2 (a 3 ) 
= u 3 (a 3 ) 



00 


TABLE II. THE RESULTS OF COMBINING SIMILAR TERMS OF THE 
LEFT AND RIGHT MEMBERS OF THE EQUATIONS IN TABLE I. 




co 


(- 1 ) 

2 


0 

h 2 (-7« 2 2 ) 

2 1 
x 0 h 2 (-a 2 2 ) + h 3 (-“3- a 2 3 ) 

3x 0 2 h 2 (-i-a 2 2 ) + 3x 0 h 3 (- “3" a; 2 3 ) + h 4 (-i-a 2 4 ) 

0 

h 2 (/3 32 a!2-^«3 2 ) 

2x 0 h 2 (/3 32 a 2 -Ao! 3 2 ) + h 3 (i3 32 a 2 2 - — a 3 3 ) 

3x 0 2 h 2 (/3 32 a 2 --|a 3 2 ) + 3x 0 h 3 (/3 32 a! 2 2 - ~s a 3 3 ) + h 4 (jS 32 a 2 3 -~a 3 4 ) 
0 

h 2 (-| a 2 2 ) 

Xoh 2 (- | a 2 V h 3 (-Tff 2 ! + a 2 2 ) 

3x oh 2 (--| a 2 2 )+ 3x(,h 3 (--|- a 2 3 + a 2 2 ) + h 4 (~ |a; 2 4 - a 2 2 + a 2 3 ) 


0 

h2 (ft 2 a 2 -~a 3 2 ) 

4 _1) 2xoh 2 (/3 32 a 2 --|a; 3 2 ) + h 3 (/S 32 a; 2 2 - Q! 3 3 - 2/3 32 a 2 + a 3 2 ) 

3x 0 2 h 2 (/3 32 a 2 --|^ 3 2 ) + 3x 0 h 3 (/3 32 a 2 2 - ya 3 3 -2(3 32 a 2 + a 3 2 ) + h 4 (/3 32 a 2 3 - ja^ + 3(5 32 a 2 - -|- 


« 3 2 -/% 2 a 2 2 + <* 3 3 ) 



APPENDIX B 


SOLUTION OF THE TWO-STEP FOURTH-ORDER EQUATIONS 


The solution of the Runge-Kutta equations for the two-step fourth-order 
process will be indicated here, since the procedure can generally be followed 
in the solution of other multi-step equations. 

We must satisfy the four equations in (12) and, in addition, the four 
applicable 12 equations in (26) , which are as follows: 


J n 


(4 ^ + <4°*) + e i2 (4 + 4 0) ) = ° ( i-B) 

[4 l) + 4 0) - 2 ] + e 12 [4 _1) (“3-1) + 4 0) «3] = 0 (2-B) 


e n 

*- \ / 

4 ^ ( e 2 i - 2ejj) + 4 ^ ( e 22 “ 2e 12 ) + 4 ^ e 2 i + 4 ^ e 22 = 0 (3-B) 


J il 


032 


4 _1) 


+ 4° A = 




where e is defined by equations (33) and (37). 

From equations (4-B) and (1-B) we easily obtain 


a; 


(- 1 ) 


+ CO 


(0) _ 


= 0 


(4-B) 


4 _1) + 4 0) = 0 


and if we combine these with the first equation in (12) , we get also 

+ cor 7 = i 


CO 


(-!) a. ,.(0) _ 


1 


( 5-B) 


(6-B) 


12 These are the equations that contain only terms in which the sums of the 
exponents in the products of the a's and /3's do not exceed three. 


19 



Using equation (5-B), both equations (2-B) and (3-B) can be reduced to 


w 2° e ll + ^3** e 12 - 0 


or 


(0) (0) o (0) ? 

u)\ ’ (2/3 Z2 a 2 ) = co 2 ’ a 2 l + co 3 a z l 


(7-B) 


Equations (5-B) and (6-B) can now be used to simplify the remaining 
three equations in (12), which become 


-co 


M) 


+ CO 


( 0 ) 


+ co> 


( 0 ) 


(0) , (0) 

co 2 o ?2 + co 3 a z = 

(0) 2 . ( 0) 2 

co 2 ’a 2 + W 3 ’ -- 


2 

_5_ 

12 

3 


( 8-B) 


This linear system in the weights can readily be solved: 


-co 


.(- 1 ) - 


6 q ? 2 <* 3 - 5 ( Q !2 + a 3 ) + _£ 

12a 2 a 3 


U) 


.( 0 ) _ 


5«3 - 4 

l2Q'2(Q'3-a' 2 ) 


- co (_1) 

- -C0 2 


JO) _ 

o; 3 - 


4 - 5 0^2 

12q' 3 (Q! 3 - a 2 ) 


(- 1 ) 

- -co 3 


(9-B) 


Using the last equation in (8-B) , equation (7-B) reduces to 
1 


032 “ 


„ ( 0 ) 

6co 3 ’ a 2 


so that equations (2-B) and (3-B) are also satisfied, provided 

/ 3 = 20's (0j 3 - Of,) 


a 2 (4 - 5a 2 ) 


( 10-B) 
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